Room-temperature electrical control of polarization and emission angle in a cavity-integrated 2D pulsed LED

Devices based on two-dimensional (2D) semiconductors hold promise for the realization of compact and versatile on-chip interconnects between electrical and optical signals. Although light emitting diodes (LEDs) are fundamental building blocks for integrated photonics, the fabrication of light sources made of bulk materials on complementary metal-oxide-semiconductor (CMOS) circuits is challenging. While LEDs based on van der Waals heterostructures have been realized, the control of the emission properties necessary for information processing remains limited. Here, we show room-temperature electrical control of the location, directionality and polarization of light emitted from a 2D LED operating at MHz frequencies. We integrate the LED in a planar cavity to couple the polariton emission angle and polarization to the in-plane exciton momentum, controlled by a lateral voltage. These findings demonstrate the potential of TMDCs as fast, compact and tunable light sources, promising for the realization of electrically driven polariton lasers.

Two-dimensional direct bandgap semiconductors are promising materials for on-chip integrated optoelectronic devices 1 . Since the demonstration of the first ultrasensitive photodetectors 2 and light emitting diodes 3,4 , there has been a strong effort to improve the quantum efficiency, speed, energy consumption and wavelength range of 2D-based devices 5,6 . In recent years, integration with CMOScompatible photonic platforms has proven the viability of the technology for the development of optical interconnects [7][8][9] . Furthermore, the integration of 2D materials with photonic cavities has enabled the demonstration of lasers 10,11 and deterministic single photon emitters 12 , main building blocks for quantum information processing.
The spin-valley properties of TMDCs have also been exploited as a new degree of freedom, with the purpose of overcoming some of the fundamental limits related to Joule heating, speed and coherence in charge-based devices. The K and K' valleys in monolayer TMDCs are coupled to circularly polarized light via the optical selection rules 13 , enabling the realization of polarization-sensitive photodetectors 14 and light-emitting diodes 15,16 . However, so far, these devices required cryogenic temperatures and complex fabrication.
Although practical applications of LEDs for information processing require fast electrical switching and low energy consumption 17 , most devices based on 2D materials operate with direct current (DC) and present very low values of radiant efficiency, due to the large contact resistance and low quantum efficiency. The control of the emission pattern is also fundamental for the coupling to photonic waveguides and detectors, as well as for illumination applications 18 . Yet, the location and directionality of light emitted from TMDCs could only be controlled by the coupling with photonic cavities and metasurfaces, where the far field emission pattern is fixed by the geometry of the structure 19 .
Here, we overcome these limitations by integrating a pulsed light emitting diode (LED) based on monolayer WSe 2 in a distributed Bragg reflector (DBR) cavity. By applying a square train of pulses to a source electrode 20,21 , we show MHz operation with the emission of ns pulses. In addition, we demonstrate electrical control of the emission angle for valley-polarized photons at room temperature. Together, our results establish a method to electrically modulate the location, directionality and polarization of 2D LEDs, with switching times ultimately limited by the pulse linewidth.

Device design and basic characterization
Our device consists of a WSe 2 monolayer encapsulated with a thick bottom h-BN and top trilayer h-BN. Monolayer graphene is used as a transparent gate electrode and top Co/Ti contacts are employed for electrical injection through the thin h-BN tunneling barrier (Fig. 1a). Figure 1b shows the optical image of the device after contact deposition (See Methods). The uniform and high-quality factor λ/2 cavity ( Supplementary Figs. 1, 2) is formed by a bottom mirror consisting of 12 pairs of Ta 2 O 5 /SiO 2 and a top mirror consisting of 9 pairs of SiO 2 /SiN (Fig. 1c). The optical properties of the mirrors are analyzed in detail in Supplementary Fig. 3. The cavity structure is designed such that monolayer WSe 2 is located at the maximum of the electric field (Supplementary Fig. 4). Figure 1d demonstrates electron and hole transport for the cavity integrated WSe 2 at room temperature (RT). The efficient electrostatic gating is confirmed by the gate-dependent PL intensity (Supplementary Figs. 5, 6). To achieve electrically driven light emission, we apply a high frequency square pulse (V sq ) with rise time τ rise < 10 ns to the source contact while we ground both the drain and the gate. The mechanism behind pulsed electroluminescence (EL) is depicted in Fig. 1e. At time T1 (Fig. 1f, inset), the monolayer is initially n-doped. When V sq switches to positive values with a rising time on the order of ns, the chemical potential of the monolayer cannot react as fast as the Fermi level of the metal, leading to a steep band bending at the interface 21 . Consequently, an electron tunneling current flows towards the contact while a hole tunneling current flows outwards from it. The spatial overlap between electron and hole wavefunctions at T2 leads to the formation of excitons and subsequent light emission. This mechanism produces the emission of short pulses with frequencies in the MHz range and duration below τ FWHM = 3 ns, as shown in Fig. 1f. The maximum operation speed can therefore reach 1=τ FWHM = 0: _ 3 GHz, limited by the device geometry and rising time of the pulse generator.

Exciton-photon coupling
Before top cavity growth, the high-quality hBN encapsulation allows us to achieve narrow excitonic emission at low temperature (Fig. 2a). At T = 5 K, we identify neutral exciton (X 0 ), charged exciton (X −/+ ), negatively charged biexciton (XX − ) and bound exciton (L) emission 20 , whereas the neutral exciton dominates the spectrum at T = 300 K. Similar results are obtained in PL measurements (Supplementary Figs. 5,6). Temperature also modifies the exciton oscillator strength and exciton lifetime, crucial for reaching the strong coupling regime 22,23 . In Fig. 2c, we plot the temperature dependence of the normalized reflectance at V g = 0 V. There is a clear correlation between reflectance and EL, with neutral and charged exciton absorption peaks appearing at low temperatures. The oscillator strength is proportional to the reflection contrast and linewidth 24 , and remains high at room temperature.
The top DBR growth results in a uniform cavity mode at the device 25 (Supplementary Fig. 7). For the bare cavity with hBN, the quality factor is Q = 296.9 ± 0.5 ( Supplementary Fig. 8). The addition of the monolayer graphene gate reduces the quality factor to Q = 179.2 ± 0.5, with a cavity linewidth γ c = 9.132 ± 0.006 meV (Fig. 2g). To quantify the coupling between WSe 2 excitons and the cavity mode, the setup is modified to image the back focal plane of the objective 26 , as indicated in Fig. 2e. The energy dispersion of the heterostructure consisting of DBR/Gr/hBN/WSe 2 /hBN/DBR is plotted in Fig. 2h. The data points are obtained from Lorentzian fits to the reflectance dispersion ( Supplementary Fig. 8). There is a clear Rabi splitting with the formation of a lower (LP) and upper (UP) polariton branches, coming from the strong coupling between the cavity mode and WSe 2 excitons. Numerical simulations considering a coupled oscillator model 27 (see Methods) allow us to reproduce the experimental results and extract a  coupling strength g = 13.4 ± 1.0 meV, with a detuning between the cavity mode and the exciton energy Δ = E c − E X = −21 meV. The proximity to metal contacts and the presence of the graphene gate reduces the coupling strength compared to previous reports 28,29 . Supplementary Fig. 11 shows further proof of the strong coupling by measuring the Zeeman splitting under external magnetic fields.

Electrical control of light emission
After characterizing the optical and electrical properties of the LED, we further apply a lateral voltage on the monolayer WSe 2 to electrically control the spatial, angular and polarization properties of emitted light at room temperature. With the application of a finite voltage on a second contact deposited on the flake, the band alignment and the total external electric field are modified. Both parameters have a direct impact on the velocity and location of carriers before recombination, which enables the electrical displacement of the polariton cloud, as shown in Fig. 3. The spatial emission profile for three representative lateral voltages V lat is shown in Fig. 3a. The cross section along the propagation direction is plotted in Fig. 3b, where d = 0 μm denotes the position of the source electrode. By performing a Gaussian fitting of the EL intensity profile, it is possible to extract the EL area and position, which is linearly dependent on V lat (Fig. 3c). We note that the conformal growth of the top mirror ensures a high cavity Q-factor even when emission occurs close to the contact edge (See Supplementary  Fig. 7). The uniformity of the cavity mode and coupling strength is evidenced by measuring the polariton emission energy along the propagation direction ( Supplementary Fig. 10).
The lateral voltage also determines the resulting emission angle of photons out of the cavity, by modifying the momentum of electrons and holes leading to the formation of polaritons. Since the movement of carriers in WSe 2 can be described by diffusive transport, their velocity is given by their mobility and the electric field: v e,h = ± μ e,h E = ± μ e,h V lat =w, with w the distance between electrodes. Since the kinetic energy of single carriers is much lower than the exciton binding energy, excitons will be generated as long as electrons and holes spatially overlap. This remains valid despite the large lateral voltages applied to the sample, since the electric field always remains lower than the dissociation field for excitons E max < E diss 30 . Considering that the scattering rate for electrons and holes is different, the exciton momentum is p X = m * e v e + m * h v h ≠0, with m * the carrier effective mass. Excitons with a finite in-plane momentum relax to the LP branch, where emission out of the cavity occurs through the photonic component, with p LP = p γ = _k == . The in-plane momentum k == is related to the photon emission angle θ by k == = 2πn c λ 0 tanθ c , with n c sin θ c = sin θ, and θ c the angle inside the cavity 31 . This implies that the polariton emission angle is linearly proportional to the applied lateral voltage to a first order approximation ( Supplementary Fig. 12). The EL dispersion as a function of θ is plotted in Fig. 4a-c for three representative values of the applied lateral voltage. The lateral voltage V 0 lat = À3 V for which EL is emitted normal to the 2D plane of the monolayer is different from zero because the dynamics of carriers is also affected by the voltage V sq applied to the source. When V lat >V 0 lat , the electrically generated excitons acquire a net momentum in the −x direction, and photons are emitted with an angle θ < 0 with the direction perpendicular to the monolayer plane. Conversely, V lat < V 0 lat leads to photons emitted with angles θ > 0. The energy and momentum conservation conditions are relaxed at room temperature (k B T~26 meV) due to phonon scattering 32 , which results in photons with a broad distribution of energy and momentum for a fixed value of the lateral voltage, as shown in Fig. 4.
Together with the photon emission angle, we can also electrically control the EL polarization at room temperature. Electrical control of EL polarization has so far only been demonstrated for electric double layer transistors based on WSe 2 . These devices, however, require an ionic liquid and can only be operated in DC at cryogenic temperatures 15 . Figure 4e-g shows the momentum and energy resolved polarization for EL at three representative values of the lateral voltage. The circular polarization can be tuned from ρ EL > 20% to ρ EL < −20%. The origin of the observed polarization comes from the anisotropic band dispersion of carriers in WSe 2 (see Methods). Under the application of a lateral electric field, the electron and hole distributions shift in momentum space. This shift induces an energy splitting between K and K' excitons, which is proportional to the in-plane momentum. Exciton relaxation to the lowest energy state in one of the valleys before scattering to the polariton can then explain the observed valley polarization (Supplementary Fig. 15). The direction of the electric field with respect to the crystal lattice is crucial for achieving a high degree of polarization. In our case, the contacts are deposited perpendicularly to the long straight edge of the WSe 2 flake, which tends to have either armchair or zigzag edges 33 . Supplementary  Fig. 13 shows the EL polarization as a function of the driving voltage amplitude V sq before the cavity growth, which supports the trigonal warping effect 15,34 as the origin of the polarization change in our device (relevant for fields on the order of 1-10 V/μm) with negligible contribution from the nonlinear valley and spin currents (with a quadratic electric field dependence). Further proof for the origin of the polarization is shown in Supplementary Fig. 14, where a second device is fabricated with contacts deposited along the armchair and zigzag directions of the monolayer. In this device, the polarization is only observable for the contact placed along the armchair direction, as expected from the trigonal warp effect. Figure 5a shows the weighted arithmetic mean of the polarization for each value of the lateral voltage (see Methods), which follows a linear behavior for moderate values of the lateral voltage. The quasi-resonant electrical injection (see Methods) together with a strong coupling to the cavity mode can explain the large degree of circular polarization shown in Fig. 5 35 . In most experiments, non-resonant optical excitations are used to create a population imbalance in the K and K' valleys. However, these high-energy excitons have a large center-of-mass wave vector. Under these conditions, there is a fast depolarization due to the long-range Coulomb exchange interaction from electrons and holes 36 . Furthermore, the L-T energy splitting (for excitons with their dipole moment parallel and perpendicular to the wave vector) increases linearly with the wave vector, leading to a fast depolarization 37   and valence band edges. In addition, the large spatial extent of the polariton wavefunction 38 reduces the intravalley scattering due to the disorder potential 39 , resulting in a longer valley pseudospin relaxation time compared to intralayer excitons 37,40 . Furthermore, the negative detuning in our cavity reduces the intervalley scattering, which only occurs through the excitonic component. To exclude the possibility that the observed polarization and directionality arises from the displacement of EL with respect to the focal point of the objective, we perform a calibration experiment ( Supplementary Fig. 16), discussed in Methods. We further demonstrate that the spatial variability of the cavity coupled exciton emission angle and polarization is negligible compared to the modulation achieved by the lateral voltage (Supplementary Fig. 17).

Discussion
The successful integration of electrically tunable light emitting diodes in planar cavities achieved in this work (Fig. 5b) constitutes a crucial step in the development of practical devices to study light-matter interaction, where the frequency, position, directionality, and polarization can be tuned electrically at room temperature. This tunability would facilitate the coupling of EL to other photonic structures such as diffraction gratings, waveguides and resonators, and can relax the alignment challenges of optical to electrical interconnects. The simplicity of the method for EL generation implies that large-area CVD grown TMDCs with metal contact arrays and top cavity growth could be used to achieve high integration densities. Furthermore, the large tunneling current per voltage pulse together with the high frequency operation could lead to a transient population inversion and electrically driven lasing. In the regime of strong-coupling, polariton lasing could be achieved at lower thresholds, with an electrically tunable chirality 41 .

Device fabrication
The substrate consists of a DBR mirror formed by 12 pairs of Ta 2 O 5 / SiO 2 layers deposited by ion beam sputtering (IBS) on top of a doped silicon wafer. Single layer graphene flakes for the bottom gate were obtained by mechanical exfoliation from graphite (NGS) on the DBR substrates and patterned to the final shape by electron beam lithography and oxygen plasma etching. Few-layer h-BN and monolayer WSe 2 (HQ Graphene) were stacked using the dry polymer-assisted transfer 42 after exfoliation on a viscoelastic stamp. The thin h-BN encapsulation and tunneling layer was transferred using a wet polymer-assisted transfer after exfoliation on a polymer double layer 43 . The thickness of the bottom and top hBN is 29 nm and 1 nm, respectively. After the stack was completed, the heterostructure was thermally annealed at 340°C under high vacuum (10 −6 mbar) for 12 h. Co (35 nm)/Ti (10 nm) contacts were defined perpendicularly to one of the WSe 2 crystal axes by electron-beam lithography and deposited by electron beam evaporation. The device was again annealed at 250°C under high vacuum (10 −6 mbar) for 12 h. Finally, the top SiO 2 spacer and nine pairs of SiO 2 /Si 3 N 4 layers were deposited by plasma enhanced chemical vapor deposition (PECVD). The second device analyzed in Supplementary Fig. 11 was fabricated following the same steps as the main device, but Ti (2 nm)/Au (80 nm) metal contacts were directly deposited on top of WSe 2 before the transfer of the hBN encapsulation layer.

Optical and electrical measurements
All measurements presented in this work were performed under vacuum at a temperature of 300 K, unless specified otherwise. Reflectance measurements were taken by focusing a fibercoupled white light source onto a 2.5 µm spot on the sample. The incident power was 90 µW. Photoluminescence was generated by excitation with a continuous-wave 647 nm laser diode focused to a diffractionlimited spot size of about 1 µm. The incident power was 200 µW. The spectral and spatial characteristics of the device emission were analysed simultaneously. The emitted light was acquired using a spectrometer (Andor Shamrock with Andor Newton CCD camera), and the laser line was removed with a long-pass 650 nm edge filter. For spatial imaging, light was collected using an Andor Ixon CCD camera after passing through a long-pass 700 nm edge filter. For polarizationresolved measurements, a rotator with a λ/2 plate was used to control the polarization incident on a λ/4 plate between the sample and the objective. A calcite beam displacer was placed before the spectrometer to separate the two measured polarizations. Angular resolved measurements were obtained by imaging the Fourier plane on the CCD camera after the spectrometer. Two convergent lenses are used to image the back focal plane of the objective, with the spectrometer slit closed to 10 µm. Electroluminescence was generated by applying a square voltage wave using an Agilent 33520 function generator. Timeresolved EL measurements were obtained with a silicon avalanche photodiode (APD) from Excelitas (SPCM-AQRH) and time correlated single photon counting (TCSPC) system (Picoharp 300).

Cavity reflectance and exciton-photon coupling
The normalized reflectance for the cavity mode is plotted in Fig. 2g. The back reflected intensity measured when exciting from the top mirror (R t ) follows the Airy distribution 44 : With 2ϕ = 2πνt RT , ν the frequency and t RT the round-trip time in the cavity. We use this equation to fit the reflectance spectra and extract the linewidth and resonant energy. The reflectance dispersion for the bare cavity and device heterostructure is plotted in Supplementary Fig. 8. The cavity mode energy is calculated as , with E ph the photon energy and n c = n SiO2 t SiO2 + n hBN t hBN t SiO2 + t hBN = 1:538 the effective refractive index based on the filling fractions of SiO 2 and hBN 45 . The reflectance dispersion in Supplementary Fig. 8 is fitted with a double Lorentzian function to extract the energies of the polariton branches. The energy dispersion is then fitted according to equation 27 : For the fitting, the exciton and cavity linewidths are fixed, and take values γ c = 9.132 ± 0.007 meV (Lorentzian fits from Fig. 2g) and γ X = 27.9 ± 0.2 meV (Gaussian fits from Supplementary Fig. 6). The coupling strength extracted from the fitting is g = 13.4 ± 1.0 meV. The Rabi splitting can then be calculated as ℏΩ = 19 ± 3 meV, with . The condition for strong coupling g > |γ X Àγ c | 2 = 9:3 ± 0:1 meV is therefore satisfied in our sample. In addition, for the two resonances to be spectrally separable, the energy splitting needs to be larger than the sum of the half-linewidths: ℏΩ > (γ X + γ c )/2 = 18.5 ± 0.1 meV. The fact that the Rabi splitting is comparable to the half-linewidths makes the observation of strong coupling difficult in our sample.
The measured splitting in reflectance is usually smaller than the intrinsic energy splitting of polariton states. In the high reflectivity limit, and assuming exciton-cavity resonance 27 : This equation can be used to recognize the underestimation of g in reflectance measurements. Considering the observed Rabi splitting of ℏΩ ref = 19 ± 3 meV and the experimental values of γ X and γ c , the coupling strength can be derived from the above equation as g i = 21.3 ± 0.6 meV. The coupling strength best describing the EL energy dispersion in Fig. 4 is g EL = 17.4 ± 0.9 meV = (g ref + g i )/2.

Pulsed electrical injection
The transient origin of current injection results in a localized emission close to the source electrode. However, as shown in Fig. 3a, the intensity maximum is not centered at the contact, due to strong lateral electric fields that deplete the electron and hole populations away from the electrode. The initial tunneling of electrons and holes can generate hot carriers with a large wave-vector and energy above the bandgap. On the other hand, electrons and holes recombine at a distance d > 500 nm from the contact (Fig. 3a). This implies that carriers need a travel time τ t = d·w/(μ e,h ·V lat ) > 5 ps to reach the recombination spot, with μ e,h taken from 46 . Therefore they relax from the excited injection condition before exciton formation occurs (scattering rate τ s ≪ τ t ) 47 .

Valley polarization of carriers under external electric fields
Lateral electric fields shift the distribution of electrons and holes in kspace, with a distribution function for carriers that can be phenomenologically described by 15 : where is the distribution of electrons and holes at the band extrema. N c,v is the normalization factor, ϵ c,v (k) the band dispersion, μ c,v (r) the chemical potential and T eff the effective carrier temperature. The relaxation time is given by τ = μ c,v m * c,v |q| , with μ c,v the carrier mobility and m * c,v the carrier effective mass. The energy dependent emission intensity resulting from electronhole recombination is given by the Fermi golden rule, which for a constant dipole moment takes the form: where δ represents the broadening. Due to the anisotropy in the valence band of TMDCs at the K and K' valleys, the resulting luminescence is circularly polarized: This mechanism is responsible for the valley polarization observed in the light emitting diode without the cavity (Supplementary Figs. 13, 14). However, this model is not valid in the strong coupling regime. The exciton-polariton valley polarization requires a more complex kinetic equation considering relaxation rates from the exciton reservoir to the polariton states, exciton-cavity detuning, coherence and valley lifetimes. A theoretical model for the case of TMDCs integrated in optical cavities has been previously developed 38 . This model can be combined with the anisotropic band dispersion of carriers to understand the observed valley polarization. In Supplementary Fig. 15 we show a simplified model in one dimension in reciprocal space to explain the effect of a lateral electric field on the exciton dispersion and polariton polarization. First, we plot the energy dispersion for single carriers at the K and K' valleys taking into account the trigonal warp effect by considering a two-band k·p model up to third order in the crystal momentum 48,49 . For simplicity, we focus on the crystal direction k x , with x being the axis along the zigzag direction of the lattice, where the electric field is applied. The exciton dispersion for both valleys at zero electric field is then simply calculated from the electron and hole energies. Even at zero field, there is a sizable splitting at large in-plane momentum. However, both K and K' excitons are degenerate in energy at the bottom of the bands. When applying an in-plane electric field, the Fermi levels of electrons and holes shift in momentum space by qτ _ E, as described before. For an applied field of E = −1 V/μm and relaxation rate 50 τ = 150 fs, the exciton dispersion becomes highly anisotropic for the K and K' valleys, as shown in Supplementary Fig. 15c. Excitons generated at k x > 0 then relax to the K' valley, which becomes the lowest energy state. These parameters represent the experimental conditions in Fig. 4a, e. The magnitude of this effect is strongly dependent on the carrier relaxation rate. Taking τ = 400 fs, the exciton energy splitting between the K and K' valleys can reach ΔE = 5 meV. This simple model can qualitatively explain the observed valley polarization, coming from exciton relaxation to one of the valleys before scattering to the polariton state.

Simulations of exciton-photon momentum and energy conservation
To theoretically compute the photon emission angle as a function of the lateral electric field, we apply momentum conservation to the recombination process (p γ = p LP = p X 0 ). Here we neglect Coulomb interactions and phonon scattering. The photon momentum can be calculated from its emission angle with respect to the normal direction to the monolayer plane θ, as described in the main text. The exciton momentum is calculated from the electron and hole velocities, p X 0 = m * e v e + m * h v h . Carrier velocities are limited by the mobility in the diffusive regime, with v e,h = μ e,h V lat /w, where μ e,h is the carrier mobility, V lat the lateral voltage and w is the width of the channel. The effective mass for electrons and holes is taken as m * h = 0:51 m e and m * e = 0:39 m e 47 . The dependence of the photon emission angle with the lateral field is calculated in Supplementary Fig. 12 for different values of the electron and hole mobilities. It should be noted that the simulations can only qualitatively describe the observed change in the photon emission angle. A more detailed description considering the carrier dynamics at the metal-semiconductor junction would be needed for a quantitative description but is outside of the scope of this work.

Reference experiment
To prove that the observed polarization and directionality is not an artifact related to the spatial displacement of EL at different lateral voltages, we perform an experiment where the objective position is scanned for fixed EL emission conditions. Supplementary Fig. 16 demonstrates that the movement of the focal point with respect to the EL location cannot explain the large change in emission angle and polarization observed in the main text.
Magnetic field dependence of electroluminescence Further proof for the presence of strong coupling in our device is given by the external response of the LP state to an external magnetic field, which is shown in Supplementary Fig. 11. The exciton and photon fractions in the LP state are given by the Hopfield coefficients 51 , which are calculated for our device in Supplementary Fig. 9. The negative detuning Δ = −21 meV implies that at θ = 0 rad emission angle, the photon fraction |C 0 | 2 is much larger than the exciton fraction |X 0 | 2 , resulting in a weak interaction of the polariton state with the magnetic field. At θ = 0 rad, the Zeeman splitting for the two circular polarization states is ΔE RÀL = g LP μ B B, with g LP = 0:41 ± 1:19 the polariton g-factor and μ B the Bohr magneton. This value is not statistically significant since the noise is higher than the energy splitting with magnetic field. At larger emission angles, the Zeeman splitting is higher due to the increase in the exciton fraction of the lower polariton. At θ = 0.15 rad, the g-factor becomes g LP = 1.06 ± 0.28. The exciton g-factor can be calculated as g X = g LP |X 0:15 | 2 = 3:53 ± 0:93, in good agreement with previous reports 52 . Exciton polaritons from the K and K' valleys have the same energy and effective mass at B = 0 T, but they preserve the circular polarization associated to the spin-orbit (SO) coupling of the bare monolayer exciton. Supplementary Fig. 11c demonstrates how EL circular polarization can be tuned with external magnetic fields.

Data analysis
Electroluminescence polarization is defined as ρ EL = ðI K À I K' Þ=ðI K + I K' Þ À ρ 0 EL , where I K (I K' ) corresponds to right (left) circularly polarized intensity and ρ 0 EL corresponds to the circular polarization obtained when V lat = 0 and V sq = 16 V.
The extracted polarization in Fig Fig. 5 corresponds to the total integrated intensity for all pixels in the detector, with the subtraction of the dark counts I T EL = ∑ðI i EL À I dark EL Þ. For magnetic-field-dependent EL measurements, the degree of circular polarization ρ EL is calculated by taking the average polarization of the N = 130 values around the peak energy, corresponding to an energy window of k B T/2, with k B the Boltzmann constant and T = 300 K. The emission energy E EL is calculated as a weighted average over same energy window, with E EL = where E i is the emission energy for a point in the spectra and I i T = I i K + I i K' the corresponding intensity. The standard errors are computed considering the error in the intensity of each pixel I K as ΔI K = 60 cts (intensity fluctuations), the error in the spectral energy E as ΔE = 0.002 meV (spectrometer resolution) and the error in the measured angle θ as Δθ = 0.03 rad (angular resolution).

Data availability
The data generated in this study have been deposited in Zenodo database at https://doi.org/10.5281/zenodo.6850668.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/.